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ABSTRACT 

We consider the extended SU(3) Nambu and Jona-Lasinio model with explicit vector couplings 
in the presence of external fields. We study the chiral anomaly in this model and its implications 
on the properties of the nucleon described as a chiral soliton of three valence quarks bounded in 
mesonic background fields. For the model to reproduce the QCD anomaly it is necessary to subtract 
suitable local and polynomial counterterms in the external and dynamical vector and axial-vector 
fields. We compute the counterterms explicitly in a vector gauge invariant regularization, and obtain 
modifications to the total effective action and vector and axial currents. We study the numerical 
influence of those counterterms in the two flavour version of model with dynamical a, tt, p, A and u> 
mesons. We find that, for time independent hedgehog configurations, the numerical effects in the 
nucleon mass, the isoscalar nucleon radius and the axial coupling constant are negligibly small. 
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1. Introduction 

The QCD chiral anomaly in its vector gauge invariant form [1] represents an important constraint for 
any effective low energy model of hadronic interactions. It is responsible for the 7r° — > 27 decay, provides 
evidence that the number of colours iV c is equal to three and is an exact one loop result not subjected to 
renormalization [2] . From the point of view of the effective Lagrangian and depending upon the degrees of 
freedom involved, there are at least two ways in which the chiral anomaly can be introduced. In a pure 
mesonic theory one can add the vector gauged Wess-Zumino term [3,4] to the action which by definition 
saturates the QCD anomaly equation. This type of construction and similar ones have been extensively used 
to calculate abnormal parity mesonic processes [5,6] or to provide short range repulsion in the construction 
of topological chiral soliton models of the Skyrme type [7]. Proceeding in this way there is in principle an 
ambiguity concerning the possible abnormal parity non-anomalous terms which do not contribute to the 
anomaly equation. One should say, however, that the ambiguity is completely removed if only pions and 
external fields are considered, but remains if other particles like e.g. heavy vector resonances are included. 
This allows to fix the corresponding coefficients in the chiral Lagrangian to fulfill some desired properties 
motivated by phenomenology and not directly derivable from the chiral anomaly [6] . 

In a relativistic chiral quark model the situation is slightly different. The non invariance of the functional 
integration measure or equivalently of the fermion determinant under chiral transformations [8,9] guarantees 
the onset of a chiral anomaly in the effective theory whose particular form depends on the symmetries 
preserved by the regulator. Actually, it is by no means clear that the resulting anomaly coincides with the 
QCD anomaly. In such a situation, the question arises whether the quark model can be redefined in such a 
way that the correct QCD anomaly can be reproduced. 

A prototype of a chiral quark model for hadronic structure is represented by the Nambu-Jona-Lasinio 
(NJL) model [10]. This model has been studied intensively in the vacuum sector, the meson sector and 
the baryon sector [11]. It is a pure quark model which incorporates explicit, spontaneous and anomalous 
chiral symmetry breaking, describes mesons as quark-antiquark pairs and baryons as solitons of three bound 
quarks or alternatively as quark-diquark bound states. However, it does not include confinement and it is 
not renormalizable requiring the use of a low energy cut-off. 

The generalization of the model to include vector and axial-vector couplings has also been studied in 
much detail either in a bosonized version [12,13,14,15,16,17] or in a pure quark language [18,19,20,21] since 
it allows to describe more mesonic phenomenology and provides a realization of vector meson dominance 
through current field identities [22]. 

Actually, we have shown in a recent paper [23] that the Nambu and Jona-Lasinio model with scalar, 
pseudoscalar, vector and axial mesons as it is most often used does not reproduce the anomaly of QCD, the 
correct result being reproduced if the vector couplings are set equal to zero. This result is perfectly consistent 
with vector gauge invariance, i.e. only the divergence of the axial current is modified with respect to the 
QCD result whereas the vector current is still conserved. Curiously, the 7r° — > 27 decay amplitude remained 
unchanged, but other anomalous amplitudes involving more than one pion, as for instance 7 — > 37r deviated 
from the current algebra result by ~ 20% for typical values of the parameters. This is an unpleasant feature 
because questions the meaning of previous calculations involving anomalous vertices, both in the meson [24] 
and the soliton sector [25,26,27,28,29,30,31] based on actions with vector mesons not containing the proper 
anomalous structure [1]. The problem arises whether the non fulfillment of the QCD chiral anomaly is a 
problem of describing vector mesons through vector couplings or whether it can be mended by addition of 
legitimate countertcrms in order to reproduce the QCD anomaly in its Bardeen form. In any case, it is 
desirable to know the possible implications not only for mesonic decays but also at the nucleon level. To 
overcome this problem Bijnens and Prades [32] have suggested, following elegant mathematical arguments 
already given by Bardeen and Zumino [33] and pursued by others (see e.g. [5] and references included in [32]), 
that indeed such counterterms exist and that they are unique if CP invariance is invoked. This represents 
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another independent non trivial solution to the anomaly equation not considered explicitly in our previous 
work [23] . The choice between these two solutions depends on our desire to include explicit vector and axial 
couplings in the starting NJL Lagrangian and at the same time being able to reproduce the QCD chiral 
anomaly. However, it should be stressed that the subtraction involves not only the external fields but also 
the dynamical vector fields. From the point of view of perturbation theory this corresponds to change an 
infinite number of diagrams in contrast to the usual external field subtraction procedure. In this respect 
this is a very particular feature of the NJL model and its singularity structure. As a consequence vector 
meson dominance in the usual sense is lost. This means that a consistent Lagrangian incorporating both 
the QCD anomaly and complete vector meson dominance cannot be constructed, at least in the vector field 
representation of vector mesons. 



In the approach of ref. [32] the calculation of the counterterms requires an explicit knowledge of the 
vector and axial currents. In addition, the formulas given by those authors are rather compact and of low 
practical utility in the present context. In the present paper we adhere to the point of view of ref. [32] 
although propose a somewhat different methodology which only relics on the knowledge of the anomaly 
itself. In fact we trace back the origin of the ambiguity to the natural regularization suggested by the 
bosonization procedure widely used for performing low energy expansions, determination of mesonic vertex 
functions and description of baryons as solitons. We derive the modifications of the effective action and 
the vector and axial currents and obtain the leading large N c corrections to the Current-Field Identities 
[22]. We also particularize the resulting formulas for the the specific two flavour case with a, n, p, A and u 
mesons and compute the numerical modification of the nucleon energy, the isoscalar nucleon radius and the 
axial coupling constant in the solitonic picture of baryons. As a byproduct we also compute the leading low 
energy chiral invariant, i.e. non anomalous, contribution to the abnormal parity action of vector mesons in 
the presence of external fields. We also study the form of the possible CP violating currents. 



The paper is organized as follows. In section 2 we review the model to fix the notation used along 
the paper. In section 3 we derive the general Ward identities from the NJL generating functional, and in 
particular, the form of the chiral anomaly which appears in the model. In section 4 the regularization in 
Minkowski space is described. In section 5 we describe the counterterms which allow to reproduce the QCD 
anomaly as well as their influence on the effective action and vector and axial currents. The constructive 
method used to derive the counterterms is described in detail in section 6. Section 7 is devoted to the study 
of the particular two flavour case with a, n, p, A and to mesons. In section 8 we present our numerical results 
for nucleon observables. Section 9 deals with the effective low energy abnormal parity but not anomalous 
action of vector mesons in the presence of external currents. Finally, in section 10 we summarize our results 
and present our conclusions. 



2. The NJL Model and Bosonization Revisited 



Our starting point is the Nambu-Jona-Lasinio Lagrangian in Minkwoski space [11] 



G f ~ X 

£njl = q(ty - M )q + -y ((q^aqf + {qKi^qf) 

G 

- — ((^t^) 2 + (qKj^q) 2 ) 

where q = {u, d,s,.. .) represents a quark spinor with N c colours and Nf flavours. The A's represent the 
Gell-Mann matrices of the U(Nf) group (see Appendix A) and M = diag(m u , m<j, m s , ...) stands for the 
current quark mass matrix. In the limiting case of vanishing current quark masses the NJL-action is invariant 
under the global U(Nf)n ® U{Nf)i, group of transformations (see e.g. Appendix A). The corresponding 
vector and axial currents are given by 

J» a {x) = \q{x)l^ a q{x)\ tf a (x) = ^q(x)j^ 5 X a q(x); (2.2) 
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respectively. We will not consider the effects of a ?7a(1) breaking term as done in [35] since they are not 
relevant for what follows. In order not to overload the paper with notation we will always work in Minkowski 
space and will never specify the Wick rotation explicitly. In fact well defined results can be obtained by 
using the customary replacement Mo — ► Mo — ie. Nevertheless all our results can be equivalently obtained 
from a Euclidean analytical continuation following the conventions of appendix B. 

The vacuum to vacuum transition amplitude in the presence of external bosonic (s,p, v, a) and fcrmionic 
(77, 77) fields of the NJL Lagrangian can be written as a path integral as 

Z[s,p, v,a,T], fj] = (0|Texp|i / d 4 x \qU + ^75 - (s + ij 5 p) Jq + fjq + qrj i|0) 

r r r r / 11 (2 ' 3) 

= / DqDqcxph / d 4 x £ N JL + qU> + ^75 - (« + *75P)J<7 + VQ + <iv\ ) 

Following the standard procedure [36] it is convenient to introduce auxiliary bosonic fields so that one gets 
the equivalent generating functional 

Z[s,p,v,a,T},fj]= J DqDqDSDPDVDA exp ji J d 4 x£ sm (x)| (2.4) 

where the semi-bosonized Lagrangian reads 

4m = £int + £-cv:t + ^cxt + + ^Mo (2-5) 

and 

Ant = q(V + 475 - (S + ij 5 P)) q = -qM iDt q 
£fL = <l(i@ + f + ^75 - (s + il5P))q = qi~D cxt q 

£L=VQ + W (2-6) 
An - tr (S 2 + P 2 ) + -±- tr (^^ + A,An 

£m = ~q~M q 

Here (S, P, V, A) are dynamical internal bosonic fields, whereas (s,p,v,a) and (r], fj) represent external bosonic 
and fcrmionic fields respectively. The bosonic fields are all of them expanded in terms of the A flavour matrices 
(See Appendix A) . Notice also that for the path integral in the bosonic fields to be well defined in Minkowski 
space we must use the prescription ^ — > — ie and — ► — ie. Finally if fermions are formally 
integrated out one obtains 

Z[s,p,7;,a,77,r7] = J DM int Det(iD) exp{-7(r7, (iD)-^)} cxpji J d 4 x£ m } (2.7) 

where the bosonic integration measure DM mt = DSDPDVDA and the usual notation 

(r?, (iD)" 1 /?) = J J d 4 xd 4 x'f](x)(x\(iD + ie)- 1 \x'}r](x') (2.8) 

have been introduced. Concerning the previous manipulations and anticipating the forthcoming sections, 
a few remarks should be made. The path integral above is a highly singular object and defines a whole 
family of effective physical theories depending on the particular prescriptions employed to give it a meaning. 
This includes in particular the regularization procedure and the symmetries respected by it. In the absence 
of a real derivation from QCD of the NJL model, one can only hope that such a freedom can be used to 
reproduce as many known features of the underlying fundamental theory as possible. This set of possible 
choices is, however, not completely arbitrary. It only reflects the ambiguity arising in the regularization of a 
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certain ultraviolet divergent Feynman diagram, namely the fact that in a loop graph the divergent piece is 
a polynomial in the external incoming momenta and also that only diagrams with less than five legs in four 
space-time dimensions are divergent. The coefficients of such a polynomial depend on the regularization. 
In practice, this means that depending on the regularization one can always subtract local and polynomial 
counterterms in the fields to the action. In this paper we are interested in the NJL model as an effective 
theory of QCD and in the possible counterterms that allow to reproduce the QCD chiral anomaly for a given 
regularization prescription. 

It should be kept in mind that, strictly speaking, the bosonization procedure is not obligatory. Nev- 
ertheless it will prove highly convenient in what follows since it reorders the diagrammatic expansion in a 
convenient way, so that N c counting rules become self-evident. More important, it is the only known device 
to describe in practice baryons as chiral solitons. Nevertheless, the bosonization may appear to impose some 
conditions on the theory, since it treats classes of infinite graphs on the same footing. This point will be 
discussed later in more detail. 

3. Ward Identities 

For the moment we will ignore the regularization and proceed formally. The resulting expressions 
will only acquire a precise meaning when the mathematical objects involved are provided with a suitable 
regularization. For simplicity we will consider the chiral limit, i.e. we set Mq = 0, since their influence is 
expected to be small in the present context. 

To obtain the form of Ward identities we decompose first the Dirac operator into an internal plus external 
field contribution both transforming homogeneously under local chiral transformations, i.e. we define 

zD = zD 0Xt - M int (3.1) 

and 

S = s + S; V=p + P; V = v + V; A = a + A (3.2) 

A local chiral rotation of the external fields induces the change D cxt — > D® xt (see Appendix A for details and 
conventions). Let us call Z 9 [s,p, v, a, 77, f)\ = Z[s 9 ,p 9 , v 9 , a 9 , if ,fj 9 ) the generating functional so transformed. 
We make a change of variables in the dynamical bosonic fields (S, P, V, A) — > (S 9 , P 9 , V 9 , A 9 ) which in the 
infinitesimal case reads 

SW = i[ev,V]+i[eA,Ai i ] 
5A»=i[e v ,A»]+i[e A ,V»] 

SS = i[e v ,S] + {e A ,P} 1 ' ' 

5P = i{e v ,P]-{e A ,S} 

Notice that the transformation is local but homogeneous for all the internal fields, i.e. no derivative terms 
appear. This is required if V,A and v,a have to transform non-homogeneously. The bosonic measure is 
invariant under these local homogeneous chiral transformations DSDPDVDA = DS 9 DP 9 DV 9 DA 9 and 
also the bosonic mass terms L 9 m — C m . This property allows to freely choose whether the bosonic fields are 
transformed or not, i.e. whether we take iD 9 = iDg Xt — M? nt or iD 9 = iD® xt — M- lnt respectively. We have 

Z 9 [s,p,w,a,77,ry] = J DqDqDM int exp{i J d i x[qiT> 9 q + C m + cpf + fj 9 q^ (3.4) 

Integrating out the fermions we get 

Z 9 [s,p,v,a,r 1 ,fj} = J DM int Det(zD 9 )cxp{-z(j7, (zD)- 1 r ? )}cxp{i J d^xC^ (3.5) 
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Note that the term (77, (zD) 1 rj) given by eq.(2.8) is invariant under the local chiral transformations specified 
above. For infinitesimal transformations we have * 

5Z[s,p,v,a,rj,f)] = J DM int s(Sp log(iD)^Det(iD) exp|-i(r?, (iD)" 1 ??)} cxpji J d 4 x£ m } (3.6) 



Decomposing the variation into its vector and axial parts as follows 

nI-i 



(Sp log(iD)) = 2i / d A x[el(x)A a v {x) + e a A {x)A a A {a 



(3.7) 



and compare with the corresponding variation before integrating out the quarks we get the following identities 



be 



J^b v c + J^b a c 



-^fabc \qhqsc + qi\isqpc 



d»J± = -A\{x)+ 



abc 



Jlb< + J%< 



-d a 



be 



qXbqSc - qi\bl5qPc 



i 

+ 2 



q\ a rj - i]\ a q 



qi\ a j5V + ¥^aj5q 



(3i 



This identity is valid under the path integral weighted with the full semi-bosonized action given by eq. 
(2.5). Moreover, it is clear that since we always obtain the Ward identities by functional differentiation of 
an effective action, the Wess-Zumino consistency conditions [3] are automatically satisfied. 



4. Effective Action and Regularization 



The determinant of the Dirac operator is an ultraviolet divergent object. Hence we have to introduce 
some regularization. In addition, since the model is non-renormalizable the corresponding cut-off has to 
remain finite, at least for the divergent pieces. This is a clear theoretical uncertainty in the model which 
has very often been neglected besides few exceptions. Nevertheless, some constraints can be imposed on the 
basis of symmetries. For the benefit of the reader we will make some digression here about the choice of 
regularization prescription at the expense of overlapping with previous works (see specially ref . [29] ) . 



At a formal level the determinant of the Dirac operator depends on the sum of the internal plus external 
fields. This circumstance is ultimately the reason for the realization of vector meson dominance through 
Current-Field Identities in the model. Naively, one might expect that after regularization the result depends 
on the sum too. This, however, docs not necessarily have to be the case [32]. For instance, one may first 
formally expand the effective action in powers of the fields and apply a regularization prescription afterwards. 
Proceeding in this way one has the freedom to regularize each vertex separately in a way that additivity is 
not fulfilled. One can add to the Lagrangian the most general local counterterm of at most mass dimension 
four which does not depend on the sum of internal plus external fields as follows 

logDet(iD) := logDet(iD) + iAT[v,a; V,A] (4.1) 



where the bar stands for a vector additive regularized fermion determinant and Ar are the counterterms. 
Such a decomposition is convenient within the chiral soliton approach to baryonic structure. In the NJL 
model, the regularized determinant represents the contribution of the polarized Dirac sea to baryonic ob- 
servables and is usually evaluated as a regularized sum of eigenvalues, thus conserving additivity. This 
corresponds to the logDct(zD) piece. In the next section we will discuss the constraints on the counterterms 
Ar. In the remainder of this section we give our precise definition of logDct(iD). 

* The total trace Sp is made out of the space-time trace, the colour trace, the Dirac spinor trace tr 7 and 
the flavour trace tr. 
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The additive contribution to the effective action can be separated into a 75-odd and 75-even part. It is 
convenient to introduce the operator 



D 5 [5, V, V, A] = 75 D[5, -V, V, -A] l5 = -D[-5, V, V, -A] 



(4.2) 



In fact, D 5 corresponds to rotate D to Euclidean space, take the hermitean conjugate and rotate back to 
Minkowski space (see Appendix B). This definition allows to separate the action into a 75-even part (normal 
pseudoparity) and a 75-odd part (abnormal pseudoparity). The former can be regularized in a chiral gauge 
invariant manner by means of the Pauli-Villars scheme [37] 



log Dot I c 



= - Splog(iD + it) + Splog(iD 5 + it) 
= i [s P log(DD 5 + it) + S P log(D 5 D + it) 
~* I Sp Cl [ lo g( DD 5 + A ? + *c) + log(D 5 D + A 2 + it) 



(4.3) 



where the Pauli-Villars regulators fulfill c = 1, A = and J2i c i = 0> J2i C « A ? = 0- For the 75-odd part we 
formally have 

^ Splog(«D + it) - Splog(iD 5 + it) 



logDet| odd 



S P log(D 2 - it) - Splog(D 2 - it) 



(4.4) 



The main difference to the 75-even part is that the sum of the eigenvalues implied in eq. (4.4) is conditionally 
convergent without the need of an explicit cut-off. The result is however unambiguous if one further imposes 
reproducing the additive anomaly GU[V, .4] (see eq. (5.4)). This is a consequence of the absence of vector- 
gauge invariant fourth order terms of abnormal pseudoparity. A practical formal expression suitable for a 
heat-kernel expansion is given by 



— — . 1 f°° da„ 

logDet odd = -7 / — Sp 
4 Jo " 



e -ia(D 2 -ie) _ e -ia(T>i-ie) 



(4.5) 



For this expression to be well defined one must take the limit e — > + at the end of the calculation. 



Finally, it should be mentioned that our separation into odd and even parts implies in itself a regu- 
larization procedure. Instead one might consider a vector additive regularization applied directly to the 
determinant of the Dirac operator, considering the operator D only. The corresponding action so regular- 
ized would differ in general from ours by vector gauge invariant, additive and chirally breaking polynomial 
counterterms whose coefficients may depend on the regularization. This feature is due to the fact that regu- 
larization and separation into odd and even parts are not commuting operations and reflects once more the 
arbitrariness in the definition of the fermion determinant. Thus, if such a kind of regularization procedure 
were used these additional counterterms should have to be subtracted. 

5. Chiral Anomaly, Counterterms and Currents in Minkowski space 



Under a local chiral rotation the variation of the Dirac determinant can be separated into two contri- 
butions. Due to the transformation properties 



tf(DD B ) = i[e v - £A75,DD 5 ] ; <5(D 5 D) = i[t v + e A75 ,D 5 D] 

5(D 2 ) = +i[t v , D 2 ] - z{{£A75, D}, D} ; <5(D 2 ) = +i[e Vl D 2 ] + i{{t Al5 , D 5 }, D 5 } 

and using the trace cyclic property, valid under regularization, we get 

S logDet I oven = 



(5.1) 



(5.2) 
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for the even piece whereas for the odd piece we obtain 



,51ogDct| odd = -i lim Sp e Al5 (e-"V 2 + e-^l) (5.3) 

Straightforward calculation using the usual heat kernel method [38] yields the following result for the in- 
finitesimal change of the regularized determinant under local chiral transformations 

tflogDet(iD) = +ij d 4 xtr[e A (x)A A (x)] = iG A [V,A] (5.4) 



where 



+ l -{T^,A a A p } + jA^F^A? + ^[D»,A v ][V a ,A ]} 



(5.5) 



with 



= M V„ - d v V^ - i[V M) V„] = i[2>„, V„] (5.6) 



Notice that the fields appearing in the expression of the anomaly are the sum of internal plus external fields. 
The QCD anomaly in its Bardeen form depends on the external fields only and corresponds to put in eq. 
(5.5) A = a and V = v or correspondingly V = and A = 0. Hence the NJL model reproduces the proper 
anomaly * if Gy — and Ar = defined by eq. (4.1). This is the solution found in our previous paper 
[23]. However, this is not the only solution. To keep the line of reasoning straight we anticipate the result to 
be derived in the next section. Another solution (unique up to CP violating terms) is given by Gy 7^ and 



J tr (6ia{F, V} + 3iF[A, V] + Aa 3 V + a 2 [A, V] + 2a{A 2 , V} + AaVaA 



24tt 2 

+ 4aV 3 + 2iDa[a, A] + iDA[a, A] + 3iDV[a, V] + 2iDV[A, V] - V A 3 - 3V 3 A^j 
where for convenience we have used notation of differential forms with the following 1-forms 

V = V^dx^-. A^A^dx 11 ] v^v^dx 11 ; a^a^dx^ (5 



(5.7) 



and 2-forms 



and 



F =^F llu dx f *dx"; DV = ^(DV) ^dx" dx u ; etc. (5.9) 



Fy. v = d^Vv - d v v,j, -i[vn, v v ) = ilD^Dv] 
(DA)^ = [D^Au] - \D v ,Ap\ 
(DV) ltu = [D ll ,V v ]-[D v ,V li ] 
(Da)^ = [Dp, a v \ - [A,,<v] 

and dx^ dx v dx a dx 13 = d A xe^ val3 . These terms satisfy that 
$(-ilogDet(iD) + Ar) =G A [v,a] 

= £ / tr { €A [ p2 + -/ + fa fl2} + A i aFa + \ {Da)2] \ } 



(5.10) 



(5.11) 



where the fermionic contribution to the effective action is regularized in a vector gauge invariant manner. 
We point out that the former eq. (5.11) corresponds to eq. (5.5) in the particular case V — and A = 0. 



* There is a sign difference with the work of Bardeen [1] due to the different convention for the Levi-Civita 
tensor £0123 (Bardeen) = +1 whereas we have e 0123 = —£0123 = +1 • 



Notice also that the counterterms do not depend on the additive combination V = v + V and A = a + A. 
Hence we will have corrections to the usual Current-Field Identities. Furthermore, the counterterms are also 
written in a manifestly vector gauge invariant fashion. 



The counterterms can be classified according to the number of external fields. For our purposes only the 
zeroth order (modification of the action Ar ) and first order (modification of the currents AJy and A J a) 
will be needed. They are 



Ar = 
Ar 1 = 



iNc 

247T 2 

iN c 
~24tt 2 



tr 



2idV[A,V] -VA 3 -W 3 A 
v(si{dA, V} - 3i{dV, A} + WAV - 2{A, V 2 }^j 
+ a(si{dV, V} + i{dA, A} + AV 3 + 2{A 2 , V}^j 
= J tr (vAJ v + clAJa^j 



+ 



(5.12) 



The first order term can be also expressed in another form if use is made of the self-consistent equations of 
motion. The total action in Minkowski space can be written as 

W = W[s + S,p + P,v + V,a + A] + &T[v,a;V,A]+W m [S,P,V,A] - (fj, (iD)-^) (5.13) 

where we have emphasized the fact that in our particular regularization the fermion determinant W = 
— ilogDct(zD) depends on the additive combinations s + S, etc. The term in the external fermionic fields 
depends also on these additive combinations. 



The former expressions can be brought into a more appealing form by considering a saddle point ap- 
proximation of the action (which in this model becomes exact in the large N c limit). In effect, if we minimize 
the total action with respect to the dynamical fields V and A (the minimization with respect to 5* and P is 
not relevant in what follows) in the presence of external fields we find that under the path integral 

and similarly for A. On the other hand, the total currents are 

SW SAT 8 . , . 

SW SAT S . -, . (r -i r\ 

= W + lnT-SV^- 1 v) (5-15) 
<5 <5 



= -J- V+ (±-±) Ar 
2G V \5v SVJ 



where in the last step we have used the equations of motion (5.14). For definiteness we will refer to them 
as the self-consistent currents. Notice that through the equations of motion the dynamical fields acquire a 
dependence on the external bosonic and fermionic fields. One can prove that the implicit dependence on the 
external fermionic fields can be neglected in the limit iV c — > oo as well as the bosonic integration. Thus we 
get in this limit 



'IT 



Jy(x)exp|i d 4 xC cxt } |0) 



(- 



2Gi 



-V + 



s_ 

5v 



_S_ 

SV 



Ar) (0|Tcxp{i J d 4 xC cxt }]\0) 



(5.16) 
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If we set the external bosonic fields equal to zero then we get that the total currents are given by 



247r 2 f V«/3 

J A l _ A , 



2G y 
24tt 2 



V 3 } + i{<9 y V a , A' 3 } + A v A a A fj + {V v V a ,A 13 } + V v A a V 13 



i{d v V a ,V } + i{d u A a , A 13 } + {V, A a A 13 } + V v V a V (i + A v V a A p 



(5.17) 



These equations represent the leading N c modifications to the usual Current-Field Identities [22] and as we 
see they are valid in the presence of external quark fields (see eq. (5.16)). In particular, they can be used 
to evaluate baryon matrix elements or form factors. It is interesting to notice that in terms of right and left 
field representation (see Appendix A) these equations can be rewritten as 



V + 



T L - 



2G V M 
2G V " 



{vrm+wivsvi 



(5.18) 



where 

v£" = d»vz-& t vz-i\y£,vz\; vr = d»vz -d*vz-i\y£,vz] (5.19) 

Let us remind that in writing the former expressions explicit use of the equations of motion has been made. It 
is interesting to note that the correction to the current field identities coincide with the difference between the 
covariant and consistent currents related to the covariant and consistent forms of the anomaly respectively 
[33]. Thus, these currents do not posses an internal anomaly. 



6. Counterterms from the Gauged Wess-Zumino action 



In this section we describe the constructive method that we have used to compute the counterterms Ar 
given by eq. (5.7). We remind that they have to be summed to the additive vector gauge invariant fermionic 
action — ilogDct(zD) so that the total sum reproduces the QCD anomaly. To do so we will proceed by 
considering the gauged Wess-Zumino term. This represents no limitation since the anomaly is saturated by 
it. In addition, it should be kept in mind that the final expression for Ar has to be a polynomial in the 
fields. 

The left-right gauged Wess-Zumino term as a functional of the additive field combinations Ar and Al 
reads (see e.g. [34] for an explicit construction) 



t^ z [u,a r ,a l ] = r wz [u] 

J tr [ArU% + z{A R , <IAr}Ur + iU^A l UArU 2 r 



48tt 2 . 

+ iU^AlUcIArUr + l -{A R U R f + A R U R - U^A l U{Ar, <1Ar} 



(6.1) 



- U^AlUArUrAr + iU^A L UA% + -{U^AlUArY 



p.c. 



where p.c. means interchanging the R and L labels and rwz[f] represents the topological Wess-Zumino 
action [4] (see Appendix D) and U a unitary flavour field WU = 1. Here the following 1-forms have been 
defined 

Ar = vr + Vr; Al=v l + V l - Ur = UUU; U L = UdU^ (6.2) 
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Under a chiral transformation the fields transform as follows 



Sv R = de R + i[e R , v R ]; Sv L = de L + i[e L , v L ] 
SV R = +i[e R , V R ]; SV l = +i[e L , V L ] 

SU = i{e L U -Ue R ); SU* = i(e R rf - rfe R ) 



(6.3) 



i.e., the dynamical fields transform homogeneously, whereas the external fields transform non-homogeneously. 
The variation gives the right-left form of the anomaly 



ST^ Z [U,A R ,A L ] = 



N c 
48tt 2 



tr 



({A R ,dA R } - iA R ^jde R 



-p.c. 



To bring this anomaly to a vector gauge invariant form we consider the polynomial action 



T^ Z [1,A R ,A L ] 



48tt 2 



tr 



-A L {A R , dA R ] + iA L A\ + -(A l Ar) 



-p.c. 



(6.4) 



(6.5) 



where 1 means the unit matrix in flavour space (U = 1 in (6.1)). The anomaly (5.4) can be reproduced by 
the effective action 

T^ z [U,A Rl A L ]^T^ z [U,A R ,A L }-T^ z [l,A R ,A L }; 6T^ z [U,A r ,A l ] = G A [V,A] (6.6) 

This action coincides with that of ref. [23] and any other action reproducing the anomaly in terms of additive 
fields (5.4) differs by chirally invariant terms from this one. Actually, the former action (6.6) is the leading 
contribution in a gradient expansion of the 75-odd part of logDct(zD) provided the non-linearly transforming 
field U is taken as the (unique) unitary part of the flavour matrices S + iP. Hence, the terms of action (6.6) 
containing the fields either U or W are not polynomial actions due to the chiral circle condition Wll = 1. 
This anomaly, however, does not coincide with the QCD anomaly [1], since it contains the dynamical vector 
and axial fields V and A. As we know the QCD anomaly depends on the external fields v and a only. To 
eliminate the internal fields dependence we propose the most general globally chiral invariant counterterm 
depending on vector and axial degrees of freedom in a way that their variation exactly cancels the dependence 
on the internal fields in the anomaly. 



T^[vr,v l ;Vr,V l ] = /" tr{ ciVrVr + C2V R dvRVR + c 3 V R v R dv R + c 4 VrV r 

+ c 5 (V R v R ) 2 + c 6 V R dv R + c 7 dV R V R v R + c 8 VrVr^ - p.c. 

Under a chiral transformation we have 

<5r£ L = J tr I [(ci + ic 2 )v 2 R V R + (ci + ic 3 )V R v 2 R - (c x + ic 2 + ic 3 )v R V R v R 
+ c 2 V R dv R + c 3 dv R V R - (c 4 + ic 6 )v R V R + (c 4 + ic 6 - ic 7 )V R v R 
+ (2c 5 + ic 7 )V R v R V R + c 7 dV R V R + c 8 Uj de R } - p.c. 

The coefficients C\ , . . . , c§ are to be fixed by imposing that 

8^ z [U,Ar,A l ]-T^[v r ,v l -,Vr,V l ])=-^ { tr [([v R ,dv R } - iv R )de R 



(6.7) 



(6.8) 



p.c. 



(6.9) 



This equation fixes all coefficients except one 

AM 

r^[v R ,v L ;V R ,V L } = - l " 



48tt 2 



tr 



3Vrv 3 r + 2iV R {dv R , v R } + cV R l (dv R - iv 2 R ) 



+ i{dV R , Vr}v r + -(V R v R ) 2 + V*v R 



(6.10) 



p.c. 
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The undetermined coefficient c vanishes if, in addition, we impose CP invariance. Indeed, under CP we have 

CP : (V R Ut, x) - -(V L r\t, -f); (V L Ut, x) - -Wit, -x) (6.11) 

where upperscript t means matrix-transposed. After some integrations by parts one can see that CP results 
in changing c — ► — c. We will take c = for the rest of this section. A brief discussion of CP violating terms 
can be found in Appendix E. 



The following action reproduces the correct anomaly (5.11) 

T^ Z [U, v, a; V, A] = r* L z [U, Ar, A L ] - T^[v R , v L ;V R , V L ] - T^ z [1,vr, v L ] (6.12) 

i.e. 6Tyf Z [U, v, a; V, A] = Ga[v,o\. However, in a vector gauge invariant calculation of the fermion de- 
terminant the chirally breaking terms, which depend on the additive combinations Ar = v R + V R and 
Al = vl + Vl, are given by eq. (6.6). Thus we define 

Ar = r* L z [i, a r , a l ] - r* L z [i, VR , v L ] - r* L K, v L] v R , v L ] (6.13) 

so that 

Y% fZ [U,v,a;V,A]=T^ z [U,A R ,A L ] + AT[v,a;V,A] (6.14) 

Thus, Ar[t>, a; V, A] are the counterterms to be added to the vector gauge invariant fermion determinant 
Det(zD), determined up to local and polynomial chiral gauge invariant combinations. These combinations 
must be formed by using the chirally covariant objects Fr, DVr and Vr in all possible combinations. For 
completeness we repeat here the argument already given in ref. [32] . If we insist on parity conservation the 
following chiral gauge invariant combinations remain 

trDV R Vl tiF R Vl tvDV R F R , tr(Wfl) 2 , trV& tr F 2 (6.15) 

together with the corresponding left-field combinations. The first term vanishes identically as can be seen 
integrating by parts. The second term is the CP violating term already considered before. The third term 
vanishes after integration by parts due to the Bianchi identity for the field strength tensor DFr = 0. The 
fourth term can be reduced integrating by parts to the second one. The fifth also vanishes identically due to 
the cyclic property of the trace. Finally, the sixth term is a topological action in the external fields which 
are assumed to have zero winding and hence vanishes. This completes the proof that the counterterms are 
uniquely given if CP invariance is invoked. 



7. Study of the two flavour case 



It is interesting to study the two flavour reduction of the model which corresponds to a NJL Lagrangian 
of the form 



G r 

C = q(i$ - M )q + — [(qq) 2 + {qTl lb qf 



2 



(Qr^ff.q) 2 + {qf^^qf 



(7.1) 



Notice that, strictly speaking, the case G2 ^ G3 is not a particular case of the model considered in section 
2, however trivial modifications can be easily implemented to the case of interest. After bosonization we get 



g.(<? + H5T + + 475) + gap - M ) q 



*o)< 



i^V 2 + t? 2 ) + \^ 2 P {r -Pn + A"- A») + \mloo^ 



(7.2) 



with G\ = g^/fj, 2 , G 2 = g 2 /{^m 2 p ) and G 3 = g^/m^. Up to the mass term, this Lagrangian is invariant under 
the SU(2)i? <g> SU(2)l <g) Us(l) group and has been studied in detail in refs. [29,30,31]. The corresponding 
baryon, vector and axial currents are 



J"{x) = q(x)j„q{x); T£{x) = ^q(x)-y „Tq{x); J*{x) = ^-q(x)-f^ 5 fq(x) 



(7.3) 
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In the notation of previous sections the reduction of the Lagrangian corresponds to take the dynamical fields 
to be ^ ^ 

v » = 2 g pP^ ■ ^+9^^ A„ = -g p A^ ■ t (7.4) 



whereas the external fields reduce to 



1 _ _ 



(7.5) 



In this model the Ward identities acquire a simple form in the case of vanishing s,p, r\ and fj external fields, 



d^jjf = 



d» Tf, = Tf, A + St A a M 



^ = Jt A ^ + # A W - ^er a %v° (d a Vp + l -v a Av - \a a A a ) 



(7.6) 



The modification of the vector gauge invariant regularized effective action can then be obtained directly from 
eqs. (5.12) giving 

Ar = ^V 

1.2- ■ .; , (7.7) 

^u»[A u A A a + 'iff A fP\ ■ A?} 
8 J 



ifau J d 4 xe^ a0 {d^ v A a ■ f + ■ A 



In this particular case the total baryon, isospin and axial currents, in the absence of external fields, become 



SW 



7 M - 
Jy ~ Sv 



SW 



Ja Sa 



(7.8) 



respectively. Straightforward calculation yields the following results for the modification of the currents as 
obtained from eq. (5.12) 



24tt 2 



9pQuj£pi/af3 



3d"(A a uj' 3 ) + 2g p uj v {A a A f ) 
3fTd° 1 ^ + 3d v p a Lo (i + g p u v {A a A A? + p° A f) 



(7.9) 



Notice that the correction to the baryon current is a total divergence, hence the baryon number normalization 
is preserved. The total self-consistent currents, i.e. the total currents evaluated by means of the equations 
of motion can be deduced from eqs. (5.17) and give 
Baryon Current 

J* = + J^Ltfe^frA" ■ f + d v F ■ P + \g p {A v A A a + 3(7 A (T) ■ A?} (7.10) 

Vector Current 



J v 



' 9 P P " + 247r> 



g ' pgiAi^fivap 



\d v A a ^ + d u uj a AP + ^gp^ifT A A f3 )} (7.11) 



AxiaJ Current 



Jt = ~— K + ^9 P 9^ v ^{d^ a f + d v fT^ + \g p u>»(A a A A? + (T A f)} (7.12) 
9 p ^~ 
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These are the corrected current-field identities in the two flavour case in leading order in N c . To conclude 
this section we notice that if the fields and vanish, there are no corrections to any single current. In 
contrast, if the field vanishes there is a correction to the baryon current not considered in previous works 
[26,27]. 

8. Numerical Results for Nucleon Observables 

Following the standard approach [11] , a baryon can be described in terms of the corresponding correlation 
function 

n B (x,x') = (0\t{b(x)B(x')}\0) (8.1) 
B(x) being a baryonic operator in terms of quark fields. We take 

B(x) = ^e Q1 - a ^1> ai - a "^ ai a 1 (x)---<z ajVcaNc (x) (8.2) 

where (ai, . . . , ajvj are colour indices, (a\, . . . , ajyj spinor-flavour indices and <J> a iv,aw c the proper com- 
pletely symmetric spinor-flavour amplitude. The exact spectral representation of the correlation function is 
obtained as usual by inserting the complete set of eigenstates of the NJL hamiltonian in eq. (8.1), namely, 

Yl B (x, x') = 9(t - t') J2(0\B(0)\B n , k)(B n , k\B(0)\0) e -^-^ k 

(8.3) 

+ (-l) N °9(t' - t) £<0|B(0)|B n) k)(B n , k\B(0)\0)e+^- x > fc 

n 

where B n , (B n ) are the baryonic (antibaryonic) states with momentum k. Further, by chosing the branch 
t > t 1 and taking the limit t — t' = T — > —zoo, the lightest baryon, and at rest, is selected in the sum, 

IL B (x,x') = (0\B{0)\B)(B\B{0)\0)e- l ™ B (8.4) 

To carry out these steps in the large N c limit, we first write the time ordered product (8.1) as a path 
integral over fermionic degrees of freedom with weight exp(i5 f NjL)- The resulting expression, in turn, can be 
obtained by appropriate functional differentiation of the generating functional Z[s,p, v, a, 77, fj] (see eq.(2.3)) 
with respect to the external quark fields rj(x) and fj(x). After bosonization and integration of the quarks 
one gets, using eq.(2.7), 

, f DM int exp(iW) n,-i iS a a'(x,x') 
TL B (x,x') = B "c$ a > "".i — — ,ih - 1 _ a ^ y ' ( 8 .5) 

where the one particle Green function S aa >(x, x') = {x\(iD + it) a a>\ x ') nas been introduced and all external 
fields s,p,v,a and 77,77 are set to zero. The limit N c — > exo drives the functional integral in the denominator 
to a saddle point bosonic configuration which describes the mean field vacuum. On the other hand, because 
there are N c factors S a , a ' in the numerator, the dominating saddle point configuration will be different from 
that of the vacuum and will depend on x and x'. Next, the limit of large evolution time T selects the minimum 
energy stationary configuration. For stationary configurations one can use the spectral representation of the 
propagator given by 

iS aa ,(x,t;x> ,t') = - t')9(e n ) - 6(t' - t)e{-e n ))^ na {x)^ na ,{x')e- K ^^ (8.6) 

n 

with ipna(x) and e n the eigenfunctions and eigenvalues of the single particle Dirac Hamiltonian H defined by 
iD — 7o(7<9 t — H) (evaluated at the stationary bosonic configuration). Choosing the t > t' branch to create 
a baryon (instead of an antibaryon) and taking the limit t — t' = T—* —zoo we get 

Il B (x,x r ) = * B (x)<i> B (x r )e^™z (8.7) 
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where the baryon wave function 



- ^•■- a "c^(g) . • - Cl^) (8.8) 

has been introduced. Here, the valence level e va i is the lowest positive eigenvalue of the Dirac Hamiltonian. 
This level is selected in eq. (8.6) by the large Euclidean time limit. Eq. (8.7) shows that spatial translational 
invariance is spontaneously broken by the mean field approximation, corresponding to the formation of a 
solitonic configuration in the presence of valence quarks. 

The baryon mass (in the large iV c limit) comes from the two contributions in the numerator and de- 
nominator, Mb — Eb — -Evac, which are given by 



E B ^E sca + AE^ + E^ + N c e v 
E mc =EZ + AE™ c + E™ c 

where the superscripts "vac" and "sol" stands for the two different stationary saddle point configurations. 
The valence term (last term) in Eb, not present in E vac , comes from the propagators S a ,a> ■ The other terms 
come from the action in eq. (5.13). The sea regularized and the pure mesonic contribution are given by 



E sc& + AE=^ logDet(iD) = ^ logDet(iD) + iAT 
E m = -j\W m 



3.10) 



respectively. At this point it is clear that any expression for the total energy involving the Dirac eigenvalues 
will necessarily depend on the additive combinations. Similarly to the energy, one body observables, like e.g. 
mean squared radius, axial coupling constant, etc. admit (at leading order in N c ) a natural decomposition in 
terms of a valence part, a pure additively regularized part and the modifications induced by the counterterms. 
This fact can be established by repeating similar steps as for the baryon mass but with non vanishing external 
bosonic fields s,p,v and a. As an example we give such a decomposition for the baryon axial coupling constant 

g N A = gt +TT + ^9a (8-11) 
i.e. in the absence of counterterms, Ag^ would vanish. 

For stationary fields the Dirac operator can be written as 

iT> = l0 (i^-H) (8.12) 

with the one particle Dirac Hamiltonian 

H = -ia • V - {V + A 0l5 ) + a ■ (V + A l5 ) + 0(S + h 5 P) (8.13) 
In the case of two flavours with hedgehog symmetry we have ( x = x/r ) 

S = g*<j{r); P = g n T ■ x<f>(r); V = g u u(r); A a = 

a . V = ?f- a ■ (x A r)p(r) ( 8 J4 ) 



a-A=^ 
2 



• r(As(r) - ^ T (r)) + (a ■ x){t ■ x)A T {r) 



If we introduce the hedgehog ansatz in the expressions for energy (7.7) and the currents (7.9) we readily find 
the modifications for the total energy and the baryon, vector and axial currents. Their detailed form can be 
seen at Appendix C. 
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Given a vector gauge invariant regularization of the fcrmion determinant one should proceed as follows. 
The total energy is obtained by adding the valence quark, the Dirac sea contribution (vector additivcly 
regularized), the bosonizing terms and the counterterms given by formula (8.9). One should then minimize 
this total energy with respect to arbitrary variations in the fields a, <f>, p, As, At and u>, and determine the 
solution of the resulting equations of motion by some iterative method until convergence is achieved. In 
the present case, such a procedure requires many iterations [29, 30, 31]. In addition, there is presently 
some concern [39,40,41] about the validity of any of the schemes proposed so far [25,28,29,30,31]. It is clear, 
however, that in cither case the proper anomalous structure has not been included since a vector additive 
regularization has always been considered, i.e. the counterterms AE are missing. Therefore, before embark- 
ing in a rather complicated calculation we make an order of magnitude estimate treating the counterterms 
as a small perturbation of the total energy. Clearly, such a procedure can be justified a posteriori if the 
modifications actually turn out to be small. 

A quantitative estimate of the corrections can be obtained by considering any of the self-consistent 
solutions available in the literature [29,30,31], since in practice they do not differ too much from each other. 
On simple dimensional grounds and using the typical sizes of the mesonic fields (<~ f v ) and their typical 
extension (~ 1 fm) for reasonable values of the parameters we estimate the counterterm to the total energy 
to be of the order of few MeV. A more systematic estimate using the self-consistent solutions of ref. [29] 
gives values for AE < 10 MeV. This small number does not stem neither from a big cancellation among 
the different terms in the expression for AE nor from cancellations between the interior and the exterior 
of the soliton. The size of the correction is to be compared to the typical soliton energies E ~ 1500 MeV. 
Thus it seems more than reasonable to treat the counterterms perturbatively. For the currents a similar 
strategy can be considered. In this case we can compare the perturbative contributions to the isoscalar 
nucleon radius and the axial coupling constant to the self-consistent ones given by formulas (C.5), (C.6), 
(C.7) and (C.10) respectively, the difference being an indication of how far are we from a self-consistent 
solution, i.e. they represent virial theorems. Again the order of magnitude estimate predicts Ar ~ 0.01 fm 
and Ag^ ~ 0.01 in any calculation scheme, much smaller than the typical values usually found (r ~ 0.8 fm 
and g^ <~ 0.5 [29,30,31]). This result is also confirmed by an accurate calculation using the solutions of ref. 
[29]. It is beyond doubt that the modification of the Nambu-Jona-Lasinio model to incorporate the proper 
QCD anomaly docs not result in large changes for the computed nucleon properties. At first sight this is 
a bit surprising since in the meson sector 20% deviations with respect to the current algebra result have 
been obtained for 7 — > tttttt [23]. On the other hand, the corrections induced by the counterterms involve 
vector mesons only (see Appendix C) whose typical extension is l/m p . Due to this the corrections in the 
nucleon currents are of order f-n/mp ~ 0.1, i.e. they are small because of the large vector meson masses. It 
is interesting to note here that similar trends have been also found in the meson sector when calculating full 
momentum dependent abnormal parity vertex functions [42]. 

9. Effective action for Vector Mesons up to fourth order in momenta 

In this section we further exploit our results to write down an effective action for vector mesons and 
external fields up to fourth order in momenta. The interesting point is that the anomaly does not imply 
the coupling strength of strong decay processes like u> — » 3ir, u> — > pir or radiative vector meson decays like 
p — > 777 etc., although the corresponding amplitudes have abnormal pseudoparity (i.e. contain an e M „ Q/ 3 
tensor). 

As we have said in section 6, up to fourth order the correct anomalous action reproducing the proper 
QCD anomaly is given by 

T^ Z [U, v, a; V, A] = T^ Z [U, A R , A L ] - T^[v, V] - I* L Z [1, v R , v L ] (9.1) 

On the other hand the action 

T^ z [U,v,a] = T^ z [U,v R ,v L ] - T^ L z [l,v R ,v L ] (9.2) 
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corresponding to Gy = and Ar = also reproduces the QCD anomaly. Therefore both actions differ by 
chirally covariant terms. A direct computation gives 



T^ z [U,v,a;V,A}=T^ z [U,v,a 
iN, 



48tt 2 



J tr {2R{F R , V R } + DV R {R, V R ] - iR{R 2 + V%)V R 



+ I(RVr) 2 + 2iUF R V R U^V L + URF R U^V L + iUV R U^V L F L (9.3) 
+ UF R RU^V L + DV R RU^V L U + iDV R [V R , U^V L U] + UR 2 V R U^V L 
+ iUV r U^V l LVl + UV R U^Vl + ^(UV R U f V L ) 2 } - p.c. 

where the following chirally covariant 1-forms 



(9.4) 



(9.5) 



R = U ] VU = UUU - iU^v L U + iv R 
L = UVU^ = UdU't - iUv 11 ^ + iv L 

and 2-forms 

Fr = dv R - iv 2 R ; F L = dv L - iv\ 
DV R = dV R - i{V R , v R }; DV L = dV L - i{V L ,v L } 

have been introduced. The transformation properties of these objects are 

5R = i[€ R ,R]; 6F R = i[e R ,F R ]; 5DV R = i[e R ,DV R ] (9.6) 

similarly for the left combinations. It is important to mention that from the point of view of an effective 
mesonic theory each term appearing in cq. (9.3) can have an arbitrary coefficient since they are separately 
chirally invariant and do not contribute to the anomaly equation. The coefficients in eq. (9.3) represent the 
particular prediction of the Nambu-Jona-Lasinio model with vector mesons so defined to satisfy the correct 
QCD anomaly. Another important point is that these terms describe low energy vector meson strong and 
radiative abnormal parity decays although they are clearly not anomalous. In particular, if the external fields 
are set equal to zero the total action is invariant under global chiral transformations, and no internal anomaly 
appears. We will not attempt to study the phcnomenological implications of the former action. Some of 
them have been considered in ref. [43] for low momenta and in ref. [42] in a full momentum dependent 
formalism, although the explicit form (9.3) has not been given. Finally as a partial check of our results, 
we consider the low energy limit of the former abnormal parity non-anomalous action. This corresponds 
to integrate vector mesons out at the mean field level or equivalently in the large N c limit. In the lowest 
relevant order we have the following equations of motion (see [35] and [29] for more details) 

V?=\{1-9 Q A )R»; V^= l -{l-g Q A )L, (9.7) 

with the quark axial coupling constant [23,35]. Moreover 

DV R = - g Q A ) (-R 2 + iF R - iU^F L u) 
DV L = - g Q A ) (-L 2 + iF L - iUFrU^ 
giving 

T^ z [U,v,a ] V,A]=Tl ?z [U,v,a] + --- (9.9) 

where the dots denote terms of order sixth at least. Therefore up to fourth order in the chiral expansion 
there arc no corrections to the effective abnormal parity action. This result is also a consequence of the 
uniqueness of the counterterms to the Wess-Zumino action (see the discussion at the end of section 6). 



(9.8) 
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10. Conclusions and Summary 

In the present work we have investigated the anomalous sector of the Nambu Jona-Lasinio model with 
vector mesons and its possible implications on the properties of the nucleon described as a system of three 
bound valence quarks in a self-consistent solitonic background of a, n, p, A and uu mesons. In most cases, 
calculations within this model have ignored the fact that it does not reproduce the correct QCD chiral 
anomaly. For the model to do so for non-vanishing vector meson fields it is necessary to modify the usual 
definition of the fermionic determinant. This can be done by subtracting suitable local and polynomial 
counterterms to the action in the vector dynamical and external fields. This represents another solution 
which was overlooked in [23], hence the conclusion there, that the NJL model with vector mesons cannot 
reproduce the correct QCD anomaly, is incorrect. These counterterms only modify the abnormal parity 
vertices at zero momentum transfer and in the chiral limit, and hence leave many meson properties such 
as meson propagators and the momentum dependence of mesonic form factors unaffected. However, there 
appear abnormal parity modifications appear in the Current-Field identities at leading order in large N c . 
In the soliton sector clear corrections appear in the nucleon mass and the axial and vector currents. For 
hedgehog profiles we have evaluated the numerical corrections to the soliton energy, the isoscalar nucleon 
radius and the axial coupling constant induced by the counterterms. We have found that they account 
for less than 1% of the total magnitude of the computed observables. As the counterterms only involve 
vector mesonic fields, and the corrections take place at zero momentum in the amplitudes, they are mainly 
sensitive to the tail of the vector meson fields only. The smallness of the counterterms in the nucleon 
might be understood due to the high vector meson masses. We conclude that the fact that the previous 
rcgularizations of the Nambu-Jona-Lasinio model do not fulfill the proper QCD chiral anomaly does not 
have practical dramatic consequences in the previous solitonic calculations. 
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Appendix A. Chiral Transformations of the Dirac operator in Minkowski space 



Wc define the effect of local chiral transformations on the Dirac operator, both in vector-axial notation 
as well as in right-left notation in Minkowski space. We follow Bjorken-Drell [44] (Itzykson-Zuber [45]) 
notation for the gamma matrices throughout. 

A.l Vector- Axial Notation 

The Dirac operator D in Minkowski space is taken to be 

iD = i$ + Y + Ji'Y5-(S + i'y 6 P) (A.l) 

with S.Py^.A^ hermitean flavour fields S — \ J2^=o * ^aS a (x), ■ ■ •, with the flavour matrices A a normalized 
as tr (A a A&) = 25 a b and \ a Xb = {d a b c + i / a fcc)A c with a, b, c = 0, . . . , iVj, — 1. Under chiral (vector and axial) 
local transformations the Dirac operator transforms as 

J) J}9 _ e +iev(x)-ie A (xh5 J^ e ~iev (x)-ie A (x)y 5 (A. 2) 

with 

€ V (x) =^4(x)A a ; e A (x) = ^e^(x)A a (A.3) 

a a 

The induced transformations on the fields are given by 

zD -» iB 9 = i$ + f 9 + A 9 j 5 - (S 9 + i l5 V 9 ) {AA) 
In the infinitesimal case this is equivalent to 



5W = [V*,e v \+i\i- A ,A'>\ 
6A» = i[e v ,Ai i ] + [Tr,e A ] 
5S = i[e v ,S] + {e A ,V} 
6V = i[e v ,V]-{e A ,S} 



(A.5) 



where the vector covariant derivative reads 

V^ = d^- iV^ (A.6) 
In addition, the dynamical and external quark fields satisfy the following transformation properties 



(A.7) 



Sq = i(e v + e A ~/ 5 )q 
5q = -iq(e v - e A j 5 ) 
Sn = i(e v - e A j 5 )n 
Sfj = -ifj(e v + e A j 5 ) 

defined to make the bilinear forms qBq and fyD -1 ?? invariant under chiral local transformations. 
A.2 Left-Right Notation 
We rewrite the Dirac operator as 

D = DflPfl + B L P L (A.8) 

with the projection operators on chirality 

Pr = 1(1+7&, Pl = \{1-1,) (A.9) 
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The right and left Dirac operators are given by 

iT> R = i@ + A R -M; iD L = i$ + A L - (A.10) 

with 



M=S + tT: M^=S-iV; 
A R = V^ + A^; Al = V»-A» 

then we have the following transformation properties under infinitesimal chiral rotations 

SM = ie L M - iMe R 

SM^ = ie R M* - iM^e L 

6A R = d»e R + i[e R ,A R ] = [D R ,e R ] 

6AZ = d»e L +i[e L ,AZ} = M,e L } 



JK> L = JI±J* (A.14) 



(A.11) 



(A.12) 



where 

e R = e v + e A ; e L = e v - e A ; V* = d„ - iA*; = - iA% (A.13) 

have been defined. Finally, the left and right currents are given by 



20 



Appendix B. Effective action, regularization and anomaly in Euclidean space 



For completeness we give below our particular conventions for Euclidean space (hatted quantities) 



— {x , x ) 



V°(x°,x l ) = iV°(x°,x i ); V l (x ,x l ) = V l (x°,x t y i W = % = (V° , V 1 ) 

d 



d Q = -id ] d l = +d t 



-If 



r 



11 dct» 



x-y = x^y^ 



-x-y 



75 = 7 7 7 7 - «7 7 7 7 = 75 



$ = j-d=-j-d = -7^ = -0; f = j-V = -fV = Y 

;0123 _ 2 —II. c 0123 _ , _ , i 

£ — £0123 — + 1 , £ — — £0123 — + 1 

ii»"*PAi i B v d a DP = -it lia , a pAi x B v C a D p 
The factor in the exponential reads 

cxpjz y d 4 xqiDq^ = exp^— J d 4 xqit)q^ 

with the Dirac operator in Euclidean space 

whose hermitean conjugate is given by 

-iD t = -i$ + f - fi% +S- i%V 



(B.l) 

(B.2) 

(B.3) 

(B.4) 
(B.5) 
(B.6) 
(B.7) 
(B.8) 
(B.9) 

(B.10) 

(B.H) 
(B.12) 



Notice that the hermitean conjugation in Euclidean space D — > corresponds to the operation D — > D 5 
in Minkowski space (see eq. (4.2)). The fermion determinant is 

Det(iD) = exp(-W) = Det(iD) = cxp(iW) (B.13) 
The vector-additively regularized real part of the fermionic contribution to the effective action reads then 



1 



ReT^ = -Sp / — $(r 



dr 



e -rDDt +e -rDtD 



(B-14) 



Here $(r) is a generalized proper-time regularization function with $(0) = <f>'(0) = and <i>(oo) = 1. The 
Pauli-Villars regularizations correspond to the choice $ = J2i °i ex p(~ T A|). The imaginary part takes the 
form 



ImW = -Sp 

2 



log(iD) - log(-iD f ) 



(B.15) 



The variation of the determinant in such regularization involves the imaginary part of the Euclidean action 
only giving 



(SlogDet(iD) = -iSImW 

iN c 



4tt 2 



J d 4 xtr{e A (x) ^[WpJV v ][iV a ,Wp] - ^A^A v AaAp (B.16) 



^A„ \iV v , iV a \A (i - l - { [*!>„ ,iT> v ],A a Ap} + ^ [i2>„ , A v ] \iV a ,Ap]]} 



where iD^ — id^ — V M has been defined. 
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Appendix C. Results for Hedgehog Profiles 



Here are the expressions for the modification of the computed quantities and the total results for self- 
consistent fields not given in the main text. 

C.l Modification of nuclcon observables 



Mean field energy 

N f°° r 11 

AE[w, p, A] = -y^9 2 p 9u J ^r 2 dru\2p(A' s - -A' T + -A T ) 

+ A{ P > + 2 -p){A s \a t ) + ^(A s + 2 -A t ) [{A s l -A T f + p 2 ] } 



(C.l) 



Isoscalar baryon density 



(C.2) 



Vcctor-isovector current 



■^J V Ya = -UakXk^9 P g^[{A S - -, \ 



s - g^i'/ 



^-2g p up{A s + \A T )} (C.3) 
r 6 j 



Axial isovector current 



(&J A Y a = ^9 P 9u J {s i a [Suj'p - MP' + ;) - l9 P u{A S - \a t ){A s + 2 -A T )] 



-3u>'p + MP' - - 2.9 P ^(p 2 - (A s - ±A T )A T ^ } 



(C.4) 



Isoscalar radius 



(Ar2>r ° = r A ^ dr{As \ at)p 



(C.5) 



Axial coupling constant 



' - Nc ~ ~ 2 i Anr 2 dr2u; 



A 9A--^9 P 9. 3 



Jo 



2(p' + y) + \9pP 2 +9p(A s - \a t ){A s + ^At) 



(C.6) 



C.2 Self-consistent nuclcon observables 



Isoscalar baryon density 



2A T\ 1 , w , 2 . 

—)- [As - -A T ){p' + -p) 
1 



- 4.9 P (^s + 3 A T ) [{As - -A T f + p 2 ] } 



(C.7) 
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Isospin current 



(J V )a 



^iak^k 



L 9 P 



24tt 2 



9 P 9w 



(A s - -A t )lu' 



/ ,i 1 m \ Atlu 1 



2 -i 



(C.8) 



Axial current 



3p 

-a^j — H -A T + ^-^g p g u ^uj p-ujp + — - -g p u[{As - -A T )A T - p Jj j 



24tt 2 



/ / P w 1 / , , 2Aj 

-Jp + U ;p' + ty + -g p Lu(As - — )(A S + — 



(C.9) 



Axial coupling constant 



9p 



V W + ^ + - A 9p{A 2 s - l -A 2 T + r/)] } 



(CIO) 
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Appendix D. Five dimensional expression for the gauged Wess-Zumino term 



In this appendix we give an alternative form for the gauged Wess-Zumino term both in its right- 
left and vector versions, T^ Z (U, vr, vl) and r^, z (C7, v, a) respectively, as five-dimensional integrals. The 
interesting feature of these formulas is that, oposite to the more conventional form, (see eq. (6.1)), the chiral 
transformation properties become more evident. 



The topological Wess-Zumino term reads [4] 



where D 5 represents a five-dimensional manifold with boundary the compactified four-dimensional space- 
time, U R = WdU is a 1-form living in D5 and U(x,s) is an interpolating unitary flavour field, WU = 1, 
satisfying U(x, 1) = U(x) and J7(x,0) a constant matrix. The following two observations will be used below. 
Since the integrand is a closed form (i.e. locally exact) the result does not depend (modulo homotopy classes) 
on the particular choice of D5. Therefore the action is well-defined modulo 2ni. In addition, in this form 
the topological Wess-Zumino term is manifestly invariant under global chiral transformations but not under 
local ones. As pointed out in rcf. [4] the traditional minimal substitution method cannot be applied in its 
standard form since the topological action is non local. In fact the existence of the anomaly prevents a full 
chiral gauging. Thus one has to resort to other methods [38]. For instance, eq. (6.1) can be obtained by 
trial and error gauging. 



in 



(D.l) 



An alternative way to obtain the gauged Wess-Zumino term is the following. We consider the minimal 
substitution rule directly in five dimensions, i.e. we make Ur — > R with R the natural five-dimensional 
extension of eq. (9.4). This requires introducing five-dimensional gauge fields as well. Clearly R(x,s) is 
chirally covariant and hence the new action is formally chirally invariant. However, since the integrand i? 5 
is not a closed 5-form the action is no longer independent of the choice of D 5 modulo 2iri. To reestablish 
the one-valucdness of the action one has to supplement the integrand with additional terms satisfying the 
following two conditions: i) the total sum has to be closed and ii) the new terms have to vanish in the 
absence of gauge fields. This uniquely determines the result. For the right-left representation we find 



T^ z [U,vr,v l } = -^ J d tr 



-^R 5 - 2RF 2 R + iR 3 F R - RFrU^F l U 



+2iv r F 2 r - v\F R - - v % 



(D.2) 



p.c. 



One should mention that the usual integration by parts is precluded in this formula since D§ has a boundary. 
One can check that the difference between the gauged and the non gauged Wess-Zumino terms is in fact the 
five-dimensional integral of an exact 5-form which is the differential of the corresponding four-dimensional 
piece in (6.1). On the other hand, the vector representation, again obtained by Bardeen's subtraction, reads 



1 wz 



[U,v,a] 



iN c 
~48tt 2 



/ 



tr 



1 

10 J 



R 5 



2RF 2 R 



iR 3 F R - RFrU^FlU 
16 



+2iaF R F L + 4iaF§, - 8a Fr - i—a 





(D.3) 



- p.c. 



Besides the fact that the last two expressions are more compact than the corresponding ones, eqs. (6.1) 
and (6.6), the chirally breaking terms are manifestly polynomial and hence the anomaly. Let us remark 
that these polynomial terms do not form an exact differential in five-dimensions and hence they cannot be 
subtracted by adding four-dimensional polynomial counterterms. In other words the anomaly cannot be 
removed. 
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Appendix E. CP-violating currents 



As already mentioned in section 6, the counterterms are unique if CP-invariance is invoked. The only 
possible CP-violating counterterms are given by 



/ 



tr 



with c an arbitrary real constant as implied by CPT. The contributions to the currents are 

TCP \ • N - 



J tr (v 



j£ p + ajf J = -ic^ J tr (dv(VA + AV) + da(V 2 + A 2 )) 



In the two flavour case we have 



B\CP 



(E.l) 



(E.2) 



(E.3) 



For hedgehog profiles all CP-violating currents vanish except the time component of the axial current 

(T\ Nc 2~ ! ( > j_P\j_( A At \IA> A T A T^\ 

( a)o = 24^2 c 9 P x a{p(p +-) + (As - —)( A s - - 3 -)) 



(E.4) 
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